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Abstract. Within the mode-coupUng theory for the evolution of structural 
relaxation in glass-forming systems, it is shown that the correlation functions 
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be presented as an asymptotic series in increasing inverse powers of the logarithm 
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model with a square-well-interaction potential. 
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1. Introduction 

Upon compression or cooling glass-forming liquids, there evolves a peculiar relaxation 
scenario called glassy dynamics. It is characterized by control-parameter sensitive 
correlation functions or spectra which arc stretched over large intervals of time t or 
frequency w, respectively. The so-called mode-coupling theory (MCT) of ideal glass 
transition has been proposed [1] as a mathematical model for glassy dynamics. The 
basic version of that theory describes the system by M correlation functions (f>q{t), q = 
1, 2, . . . , M, which have the meaning of canonically defined auto-correlators of density 
fluctuations for wave vector moduli q chosen from a grid of M values. The theory 
deals with a closed set of coupled nonlinear equations of motion for the 4)q{t). The 
coupling coefficients in these equations are determined by the equilibrium structure 
factors, which are assumed to be known smooth functions of control parameters like 
temperature T or density p. The solution of the MCT equations describe a transition 
from an ergodic liquid state to a non-ergodic glass state if the control parameters 
pass critical values Tc or pc, respectively. This transition is accompanied with the 
appearance of a dynamical scenario, whose qualitative features can be understood by 
asymptotic solution of the equations for long times and control parameters close to the 
critical values. The asymptotic formulas establish the universal features of the glassy 
dynamics described by MCT. On the basis of this understanding, one can construct 
schematic models. These are based on equations of motion which have the same 
general form as those derived within the microscopic theory of liquids, but use the 
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number M of correlators and the coupling coefficients as model parameters. Thereby, 
one gets simplified models whose results can be used for the analysis of data [2] . 

The ideal liquid-glass transition described by MCT is a fold bifurcation exhibited 
by the equations of motion. One can show, that all singularities that arc generically 
possible, are of the cuspoid type [3,4]. Using Arnol'd's notation [5], an Ai is a 
bifurcation which is equivalent to that for the roots of a real polynomial of degree 
I. Simple schematic models using only a single correlator exhibit besides the fold 
singularity also the cusp singularity A3 and the swallowtail singularity A4 [2]. The 
bifurcation dynamics near a higher-order glass-transition singularity Ai, I ^ 3, is 
utterly different from the one near the liquid-glass-transition singularity of type A2. A 
major new feature is the appearance of logarithmic decay yielding to a much stronger 
stretching than known for the A2-scenario [6] . The result for the correlators of M = 1 
models has been worked out in a certain leading-order multiple-scaling-law limit [7]. 
These formulas have been used to fit dielectric-loss data of glassy polymer melts [8- 
11], thereby providing some hint that the MCT for higher-order singularities might 
be of relevance for understanding glassy dynamics. 

Let us consider a system of spherical particles interacting via a steep repulsion 
potential characterized by a diameter parameter d, which is complemented by an 
attractive potential. The latter shall be characterized by an extension length A and an 
attraction-potential depth uq. Such system is specified conveniently by three control 
parameters: the packing fraction ip = p-KcP /Q, the dimensionless attraction strength 
r = uo/{kBT) or the dimensionless effective temperature 9 = l/F, and the relative 
attraction width 5 = A/d. If (S is sufficiently large, this potential is a caricature of 
a van-der-Waals interaction. One gets a decreasing r'^-versus-(p'^ line of liquid-glass 
transitions in the T-ip plane of the thermodynamic states similar to what was first 
calculated within MCT for Lennard- Jones systems [12]. For F = 0, the mentioned 
line terminates at the critical packing fraction (/'hss ^'-'^ vitrification of the hard- 
sphere system. The decrease of cp'^ with increasing F'^ expresses the intuitive fact that 
cooling stabilizes the glass state. If i5 is sufficiently small, however, there appear two 
new phenomena. First, for small F, the F'^-versus-(/j'^ line increases. Cooling stabilizes 
the liquid because bond formation creates inhomogeneities which favor fluidity. The 
new liquid state for ip > ip^gg exhibits a reentry phenomenon. Classification occurs 
not only by increasing F, i.e. by cooling, but also by decreasing F, i.e. by heating. 
Second, the F'^- versus- (^'^ line can consist of two branches that form a corner. The low- 
F-branch is terminated by the the high-F-branch. The latter continues into the glass 
state as glass-glass-transition line, which has an As-singularity as endpoint. These two 
phenomena have been found by using Baxter's model for the structure factor as input 
of the MCT equations [13, 14]. In these calculations, the wave- vector cutoff (?max used 
in the MCT model defines the range parameter S = tt / {qnia,xd) ■ The generic possibility 
for a transition from small-5 states with Aa-singularity to large-(5 states without A3- 
singularity is the appearance of an A4-singularity for some critical value S* . Since the 
Ai bifurcations deal with topological singularities, the indicated scenarios are robust, 
i.e., they occur for all potentials of the kind specified above. The A4-singularity was 
identified first for the square- well system (SWS), i.e., for a system where a hard- 
core repulsion is complemented by a shell of constant attraction strength uq. Here, 
6* « 0.04 was calculated [15] . The neighborhood of the was analyzed for some other 
potentials with the conclusion that there are no qualitative differences between results 
referring to different shapes of the potential or to different approximation schemes for 
the structure factor [16]. 
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The above described systems with short-ranged attraction can be prepared as 
colloidal suspensions. The liquid-glass-transition lines can be identified by analyzing 
the nucleation processes. Light-scattering experiments can provide the density- 
correlation functions 4>q{t). Such studies have identified the existence of liquid states 
for Lp > f^gg and the reentry phenomenon [17-19]. Molecular-dynamics simulation 
studies can determine the mean-squared displacement and the diffusivitics with good 
accuracy. These quantities exhibit drastic precursors of the liquid-glass transition. 
Several simulation studies [18, 20-22] have confirmed the predictions on the reentry 
phenomenon. Near the corner formed by large-F and the small-F transition lines, there 
should occur an almost logarithmic decay of the density correlations 4>q(t), which 
is followed by a von-Schweidler-law decay as beginning of an a-relaxation process 
[13,15]. Such scenario was first reported for micellar solutions [23]. This signature 
of the dynamics for ip > (Pngg states was detected also for colloidal suspensions with 
depletion attraction [19]. 

In order to identify a higher-order singularity in data from experiment or from 
molecular-dynamics simulation, one has to identify the features of the correlators (jjqit) 
which are characteristic for these singularities. The general theory of the logarithmic 
decay laws caused by an Ai for 1^3, has been developed and the relevant general 
scenarios have been illustrated for schematic models [24, 25]. The specific implications 
of the general theory for the SWS, in particular the change of the features with 
changes of the wave number and the peculiarities expected for the mean-squared 
displacement, have been worked out as well [26,27]. Simulation data for the tagged- 
particle-density correlators as function of the wave vector q [28] provide a first hint that 
the predicted logarithmic decay processes for ip > (^^hss states near an Aa-singularity 
are present. Major progress was reported recently for simulation studies for two 
states of a binary SWS [29]. The logarithmic decay and its expected deformation 
with wave-number changes has been detected convincingly. The identified amplitudes 
agree semi-quantitatively with the calculated ones [26]. In addition, the mean-squared 
displacement exhibits the expected control-parameter dependent power-law behavior. 
These findings provide very strong arguments for the existence of a higher-order glass- 
transition singularity. One concludes that the cited MCT results on simple systems 
with short-ranged attraction reproduce some subtle features of glassy dynamics so 
that further studies of these systems within that theory seem worthwhile. 

If one shifts control parameters towards the ones specifying a higher-order 
singularity, the time interval for logarithmic decay expands. But, simultaneously, 
also the beginning of the time interval shifts to larger values. Consequently, there 
opens a time interval of increasing length between the end of the transient and the 
beginning of the logarithmic decay. Within this interval, the correlators are close to 
the critical ones (f>q{t), i.e., to the correlators calculated for control parameters at the 
singularity. It should be expected that these critical correlators will be detected in 
future data from experiments and from simulation studies. It was shown for one- 
component schematic models that the critical correlators approach their long-time 
limit proportional to 1/ ln'"(t/to), where to = 2/(/ — 2) for an Ai [7]. In the following, 
these results shall be extended in two directions. First, the critical correlators shall be 
expanded in an asymptotic series so that an estimate of the range of validity of various 
asymptotic formulas is possible. Second, the (pgit) shall be calculated for the general 
theory so that a discussion of the g-dependent corrections of the leading asymptotic 
formulas is possible for an A3- and an A4-singularity. 

The paper is organized as follows. In Sec. 2, the general starting equations for an 
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asymptotic discussion of critical relaxations arc compiled. Then, in Sec. 3 and Sec. 4, 
the asymptotic expansion is carried out for one-component models for states near an 
A3- and A4-singularity, respectively. The results will be demonstrated quantitatively 
for schematic models. Section 5 shows how the theory for (f>q{t) for an A3 can be 
reduced to the theory of one-component models, The results are demonstrated for 
a two-component schematic model and for the SWS. The analogous results for an 
A4-singularity are presented in Sec. 6. A summary is formulated in Sec. 7. 

2. General equations 

2.1. Equations for structural relaxation at glass-transition singularities 

Within the basic version of the mode-coupling theory for the evolution of glassy 
relaxation (MCT), the system's dynamics is described by M correlators 4)q{t), q = 
1, . . . ,M. The theory uses the exact Zwanzig-Mori equations of motion. These are 
specified by M characteristic frequencies fig > and M fluctuating-force kernels 
Mq(t). The latter are decomposed into regular terms M^'^^{t) describing normal-liquid 
effects and in mode-coupling kernels ruq (t) . The essential step in the derivation of the 
theory is the application of Kawasaki's factorization approximation to express the 
kernels mq{t) as absolutely monotone functions Tq of the correlators. These functions 
depend smoothly on a vector V of control parameters like density and temperature. 

mq{t)=Tq[Y,c^k{i)]. (1) 

Vector V specifies the equilibrium structure functions of the system. Using Laplace 
transforms of functions of time, say F{t)^ to functions defined in the upper plane of 
complex frequencies z, F{z) = i Jq dt exp(izi) the equations of motion read 

(j)q{z) = —l/{z — fl'^[z-\-M^^^{z)-\-{l,'^mq{z)]} [2]. Glassy dynamics is characterized by 
long-time decay processes that lead to large small-frequency contributions to mq{z). 
These small-z contributions to mq{z) dominate over z + M^°s(z). Therefore, glassy 
dynamics is described by the simplified equation 4>q{z) = —l/[z~ \/mq{z)\ [30]. 
Equivalently, there holds 4'q{z)/[\ — z(j}q{z)\ = mq{z). It will be more convenient to 
modify the Laplace transform to another invertible mapping S from the time domain 
to the domain of complex frequencies according to 

/•oo 

S[F{t)]{z)^{~iz) I dt cxp{izt) F{t) . (2) 
"'0 

Using this notation, the MCT equations for the small-frequency dynamics read 

s[Mtm/{i - s[Mm^)} = s [Tq[v, Mm (^) . (3) 

Since !Fq is determined completely by the equilibrium structure functions, the 
dynamics obtained from Eqs. (1) and (3) is referred to as structural relaxation. These 
equations are scale invariant: if <f>q(t) is a solution, the same is true for 4>q{t) = (j}q{xt) 
for all a; > 0. The scale for the dynamics is determined by the transient motion. The 
latter is governed by fi^ and Mg®8(t). Since these quantities do not enter Eq. (3), the 
solutions of Eqs. (1) and (3) are fixed only up to an overall time scale [30]. In the 
following, this time scale will be denoted by t^. 

A glass state is characterized by non- vanishing long-time limits of the correlators: 
limt^oo (t>q{t) = fq, < fq < 1- Equivalcutly, one gets lim^^^o S[(l)q{t)]{z) = fq. Hence, 
the zero-frequency hmit of Eq. (3) yields /g/(l — fq) = ^^[V, /fe], g = 1, 2, . . . , AI. 
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This is a set of M implicit equations to be obeyed by the M numbers fq [1]. If the 
Jacobian of these equations is invertible, the solutions vary smoothly with changes of 
V. If the Jacobian is singular for some state with fq = fq, fq exhibits a singularity 
as function of V for V tending towards V"^. Therefore, such state is called glass- 
transition singularity. The solution for the correlators for V = is referred to as 
critical correlator </'g(i). Let us introduce a functions 0<j(t) by 

0^W=/,^ + (l-/,^)<^,W (4) 

obeying hmt^oo 4>q{t) = 0. In the following, and S[(j)q{t)] (z) shall be used as small 
quantities for an asymptotic expansion of (/'^(t) for large times and small frequencies. 
Introducing the coefficients 

4tU. = ^ (1 - ^'){5"^Jv^ f^yofk, ■ ■ ■ dfkJii - J • ■ • (1 - /L ) , (5) 

Eq. (3) can be rewritten as the set of equations of motion for the (l>q{t) [2]: 

[Sqk~A['^']S[Mtm = M^)- (6) 
Here, Jq{z) = X]ri>2 Jq^^'i-^) with the nth order expansion term given by 

4"^ (^) = 4t.kM^. (0 • • • k„ {tm - simn^) ■ (7) 

In Eq. (6) and (7) and in all the following equations, summation over pairs of equal 
labels k is implied. The M x M matrix [6qk — A'qk'^] is the Jacobian mentioned above. 

(l)c 

Therefore, a singularity is characterized by matrix A^j, to have an eigenvalue unity. 
It is a subtle property of the MCT equations that this eigenvalue is non-degenerate 

flic 

and that all other eigenvalues of A^^. have a modulus smaller than unity. The left 
and right eigenvectors shall be denoted by a* and a^, q = 1, . . . , M, respectively, 

alAif = a*q, 4tak^aq. (8) 

Generically, one can require aq ^ 0, a* ^ for q = 1, . . . , M. To fix the eigenvectors 
uniquely, two normalization conditions can be imposed: o,qCiq ~ 1, X]<j '^g^g'^g = 1 
[2,3]. 

Because of the non-degeneracy mentioned, the singularity is topologically 
equivalent to that of the zeros of a real polynomial of degree 1,1 — 2,3,.... It is 
a bifurcation of type Ai [5]. The singularity can be characterized by a sequence of 
real coefficients /i2, /is, .... An Ai is specified by /x„ = for n < I and /i; ^ 0. The 
simplest of these numbers reads 

M2 = 1 - X! «9^9fcifc2«'=i«fc2 ■ (9) 

For an A2-glass-transition singularity, fi2 determines the so-called critical exponent 
a, < a < 1/2. In this case, the critical correlator can be asymptotically expanded as 
a power series: (l)q{t) = aq{to/t)°' -\- a'q{to/t)^'^ -I- • • • . If the A2 singularity approaches 
a higher-order singularity Ai, I > 3, the exponent a approaches zero and the cited 
asymptotic expansion breaks down [2] . It is the goal of this paper to derive a long- 
time expansion of the critical correlator at A3- and A4-singularities. Equivalently, it 
is the aim to solve asymptotically Eqs. (6) and (7) for (j)q{t) for states with 

M2 = , fisy^O (10a) 
for an Aa-singularity denoted by V = V° and 

M2=M3=0, //4 7^0 (10&) 

for an y44-singularity denoted by V = V* . 
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2.2. Expansions of slowly-varying functions 

The derivations in this paper shall be based on an extension of the Tauberian theorem 
for slowly- varying functions, which has been introduced in Ref. [7]. A function C{t) 
is called of slow variation for long times if limT-^oo C{tT) /C{T) = 1 for all t > 0. 
This is equivalent to 7(2) = S[C{t)]{z) being slowly varying for small frequencies: 
limT^oo 7(2/r)/7(i/T) = 1. In addition, the Tauberian theorem states, that 7(2) 
is asymptotically equal to G{i/z) : Ynnz-,o^{z)/G{i/ z) = 1 [31]. Typical examples 
for functions of slow variation are pm(ln(lni))/ ln'™(t), where m — 1, 2, . . . and pm 
denotes some polynomial. The critical correlator shall be expressed as sum of 

such functions. Let us introduce the notations 

G{t) ^ g{x) , x = lii{t/to), y = ln(i/zto) , (11a) 

la 

gra{x) =Pm{x)/x™ , Pm{x) = ^ Cj„,ix' . (11&) 

Z=0 

gm.jri{x) is asymptotically negligible compared to gm{x): lim2,_,oo gm+i{x) I gm{x) = 0. 
For later convenience, let us write f{x) = 0(1/0;™) if /(a;)a;™ is bounded for large x by 
some polynomial of Ina;. Denoting derivatives by d"-g{x)/dx"' = ^'■"•'(a;), n = 0, 1, . . ., 
one finds 

gL")(x) = 0(l/x™+"). (12) 

Equation (2) can be rewritten as S[G{t)]{z) = cxp{— u)g{y In u) du. Formal 
expansion in powers of In u leads to 

5[G(i)](^) = £^r„g(")(y). (13) 

Here, r„ = r'")(l) denotes the nth derivative of the gamma function at unity. One 
gets Fq ~ 1, — Fi = 7 is Euler's constant, and F„ for n > 2 can be expressed in terms 
of 7 and Riemann's zeta-function values C(^) with K = 2, . . . , n [32]. For example, 
F2 — F^ = C(2)- Using Eq. (13) with G{t) = gm{x), one gets an asymptotic expansion 
in terms of increasing order 0(l/y™+"). The leading n = contribution is gm{y)', 
and this is the result of the Tauberian theorem [31]. The terms for n > 1 provide 
systematic improvements for large y, i.e. for large times or small frequencies [7]. 
If one uses Eq. (13) for G{t) ~ G{t)F{t), one gets the asymptotic expansion 

s[Git)Fmz)-s[Gmz)s[Fmz) = Y.Y. L_ g^"'"H^/)/^"H^)-(i4) 

n=2 m=l ^ ' 

Let us use G{t) = g^i {x) and F{t) = g^^ (x). The Tauberian theorem imphes that the 
leading contribution to S[Gmi{t)Gm2it)]{z) cancels against the leading contribution 
to S[Gmi{t)]{z)S[Gm2{t)]{z)- The tricks underlying the asymptotic solution of the 
MCT equations at a higher-order singularity are based on the observation that also 
the leading corrections to the Tauberian theorem cancel [7]: 

S[g,,,At)gmAm^) - S[gmAm^)S[gm2{m^) = 0(1/2/"^+'"^+') . (15) 

The difference between the two terms on the left-hand side is two orders smaller for 
vanishing frequencies than each of the terms separately. 
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3. Critical correlators for one-component models at an Aa-singularity 

3.1. The leading contribution 

It will be shown in Sec. 5 how one can reduce the problem of solving Eqs. (6) and (7) 
for a general number M of the correlators to the special problem of solving for M = 1 
models. Therefore, the problem shall be discussed first for one-component models. 
For this case, one can drop the indices in all formulas of Sec. 2.1. There is only one 
correlator 4>'^{t), one long-time limit for the critical point V'^, and one function (j>{t) 
determining the critical correlator as (j)''{t) ~ + {1 — f'^)(t>{t). The Jacobian matrix 
agrees with its eigenvalue, and this is zero. Hence, Eqs. (6) and (7) can be noted as 

X(z)=0, (16a) 

OO 

= ^if„(z). (16&) 

Here, Kn{z) is the expansion term of order 0". Let us introduce the abbreviation 

M^) = SW\t)]{z) ^ SmTiz) , (17) 

and denote its inverse transform by ij^nit), i.e., iS['0„(i)](z) = ijjniz). Remembering 
that for M = 1 models there holds ^„ = 1 — A*^"^^, one gets K„{z) = ipn{z) — 
/x„5[(/)"(t)](z) [24]. Specializing to the As-singularity as noted in Eq. (10a), the 
equation of motion (16a) is defined by 

K{z)= V2(^)-/^35[03(i)](z) 

+ KV3(^)-/^45[04(i)](z) (18) 

+K'iz) . 

Here, K'{z) = K'ip4{z) — ^5S[(j)^ {t)]{z) + . . .. The numbers k and k' have been 
introduced for later convenience. For the AI = 1 models under consideration, one 
has to substitute k = k' = 1. 

Let us examine whether one can solve the equations with the Ansatz 0(t) = 
grnix) = Cralx^. From Eq. (13) one gets Si^'^ifyXiz) = {cm/y"'? + 0(l/y3™+i). 
Using Eq. (14) with G{t) = F{t) = g„,{x), one obtains ^2(2) = C,{2){mc„, / y™+^f + 
(J[\ly'^i^+^y Choosing m — 2, both terms in the first line of Eq. (18) are of the 
same order 1/y^. They cancel in this leading order if /^acf = 4C(2)c2. From Eqs. (13) 
and (15) one infers that the terms in the second line of Eq. (18) are of order 1/y* and 
K' ~ 0(1/?/^"). One concludes that the leading asymptotic behavior of the critical 
correlator for large times is described by (j}{t) = 32(2^), where 

g2{x)=C2/x\ C2=4C(2)^3. (19) 

3.2. The leading correction 

Let us split the function (j}{t) into its leading term and a correction g{x): 

4>{t) = g2{x) + g{x) . (20) 

Substitution of this formula into the first line of Eq. (18), one gets expressions up 
to third order in g. The term independent of g is S[g2{x)]{z) — S[g2{x)Y'{z) — 
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/Z35[g|(x)](z), and it shall be denoted by [(4C(2))Vm3]-F(2/). Equations (13) and (14) 
are used to derive the asymptotic series 



t-< M3y^+" 130 {n~2)l 



n— 3 

- ^ (n - m + l)(m + l)(r„ - r„_™r™) I . (21a) 

m=l J 

The term linear in g is 2{5[52(a;).g(x)](z)— 5[(72(a^)](z)5[g(x)](z)} — 3/X35[g|(a;)g(x)](z). 
It shaU be denoted by [(4C(2))2/^3] [Pg(y) + V'g{y)]. Here, the differential operator 
D yields the leading contribution 

■Dg{y) - [y ■ dg{y)/dy + 3g{y)]/y^ . {21b) 
The correction V is expanded with the aid of Eqs. (13) and (14): 

oc n— 1 

V'g{y) = [1/2C(2)]^ ^(-1)"-" {[g^''"Hy)/y-+'~''^m\]{r^ - T^.^T^) 

n—3 m— 1 

+C(2)r„_2[ff('"-'H2/)/2;"+'-™(m - l)l]{n -m + l)in - m){n - m + 1)} .(21c) 

With these notations, the equation of motion for g{y) is reformulated as a linear 
differential equation with some inhomogeneity I{y): 

V~g{y) = I{y) , (22a) 

I{y) = F{y)+V'~g{y) 

+ S[g\x)\{z) - S[g{x)?{z) - ii,,S[g2{x)f{x)]{z) 

- ^i^S[f{x)]{z) + nij^iz) - ^JiiS[^\t)]{z) + K'{z) . (22&) 

It might be adequate to emphasize, that Eqs. (19)-(226) formulate an exact rewriting 
of Eq. (3) for M = 1 models. 

The iterative solution of Eq. (22a ) for g{x) is based on the observation, that one 
gets for functions gm{y) from Eq. (11&): 

T^gmiv) - KAy) + (3 - ™)p™]/y'"+^ . (23) 

If one tries with g{x) = gaix), one finds on the one hand 'Dg„i{y) = p'^{y)/y'^. 
On the other hand, one verifies, that all terms on the right hand side of Eq. (22&) 
are 0{l/y^) except for the n = A contribution to F{y). One checks that F{y) ~ 
24C(3)/(/i3y^) + 0{l/y^). Hence, the leading order solution for g reads 

33 (x) = C3 ln(a;)/x3 , C3 = 24C(3)//i3 . (24) 

Combining this finding with Eqs. (19) and (20) and eliminating all the abbreviations, 
one reproduces a result of Ref . [7] : 

0°(t) = /° + (1 - /°)[c2/ln2(t/to)] {1 + [6C(3)/C(2)]lnln(t/to)/ln(tAo)} • (25) 
This formula describes the critical correlator up to errors of the order 1/ ln^(t/to)- 
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3.3. Higher- order contributions 

The equation for g{y) allows for an iterative solution so that the iteration step with 
number m reads g = (?3 + .94 + • • • + .9™- Here the numerator polynomial in Eq. (11 &) 
is of degree not larger than (m — 2), i.e., 

m-2 

= 5^ c„,,zln'(a;)/x". (26) 

Suppose, the procedure had been carried out up to step m — 1, m = 4, 5, . . .. Then 
'^9{y) = T^9m{v) + 0{l/y'^^^). By construction, all terms up to order (m + 3) cancel 
against the one appearing in I{y). One checks, that the leading contribution to /(y) 
reads p(lny)/2/'""'"*, where the degree of the polynomial p does not exceed m — 3. 
Hence, Eq. (22a) is equivalent to the linear differential equation p'^^ + (3 — m)pjn — p. 
It is readily solved by Eq. (26), provided the coefficients Cm,i are chosen properly. 

In order to determine (74 and g^, one can drop K'{z) in Eq. (22&). The coefficients 
Cm,i are given by fi^, k, and fi^ as follows 

C4,o = 792C(3)V(^V3) + [4^4/(9/^3) - 4K/(3Ai3) - 7/6] TT^fi^ , (27a) 
C44 = - 432C(3) V(7rV3) , {27b) 
C4,2 = 648C(3)V(^V3) , (27c) 



cg^o = C(3)7r^[400K/X3 + 1551/4 - 160^4/ il5i4) 

- [39744 Ci^f/TT^ + 528C(5)]/Ai3 , (28a) 

C54 = 64800 C(3) V(7rV3) - 4C(3)7r2(21/i2 - 24^1/13 + 8/^4)/^! , (286) 

C5,2 = - 27216 C(3)V(^V3), (28c) 

C5,3 = 15552 C(3)V(7rV3), (28rf) 



The coefficients for gg and have also been determined . The only new model 
parameters entering the coefficients arc fi^ and k' [33] . 

3.4. Discussion 

The preceding results shall be demonstrated quantitatively for the simplest model 
exhibiting a generic j43-glass-transition singularity. This model was derived for a 
spin-glass system and it is defined by the mode-coupling function [34] 

m{t) ^vi<l){t) +V3(f>^{t) . (29) 

Here and in the following models, we use a Brownian short-time dynamics as specified 
by the equation of motion 

Tdt(j){t) + 4>{t)+ I dt'm(t-~t')d[(l){t')^Q, (30) 
Jo 

to be solved with the initial condition 0(t — > 0) = 1. The short-time asymptote is 
4){t) — 1 = — (i/r) + ©((t/r)^). The singularity is obtained for the coupling constants 
1"° ~ 9/8 and — 27/8. The critical long-time limit of the correlator is /° ~ 1/3 
[6,24]. The other parameters entering the coefficients via Eqs. (3.3) and (3.3) are 
/X3 = 1/3 and ^4 = /j.^ = k = k' = 1. Thus, all expansion formulas are specified, except 
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for the time scale tg. To case reference to various degrees of asymptotic expansions, 
let us introduce the abbreviation for the nth order approximation 

n 

<t>°{t)n = /° + (1 - n Gnit) , Gn{t) = ^ gMt/to)) ■ (31) 

m=2 




-1 0123456789 10 

logio(t/T) 

Figure 1. Critical decay at the As-singularity of the model defined by Eqs. (29) 
and (30). The full line shows the solution for cj>°{t). The lines labeled Gn, 
for n = 2, 3, 5, 7, show the approximations from Eq. (31) with the time scale 
to/r = 1.6 ■ 10-*. The time where G3, G5, and G7 deviate by 2% from (j>°{t) is 
marked by a square (□), a triangle (A), and a diamond (o), respectively. 

Figure 1 exhibits 4>°{t) as obtained from Eqs. (29) and (30) for the state V = V°. 
The approach to the critical plateau f° is significantly slower than the one for a 
typical A2-singularity. In the latter case, the decay comes close to the plateau within 
a few decades of increase of the time when a deviation of 5% is used as a measure. 
Such a criterion is not met by the decay in Fig. 1 for the entire window in time 
shown. For t = 10^^, the critical correlator (t>°{t) is still 5.5% above f°. To apply the 
asymptotic approximation, one has to match the time scale to at large times. A reliable 
determination of is not possible when using only G2 {x) or G3 {x) . Using Gy (x) and 
extending the numerical solution to t/r = 10'^®, it is possible to fix to/r = 1.6 • 10"''. 
Notice, that to is several orders of magnitude smaller than the time scale t for the 
transient dynamics. With this value for t^, the successive asymptotic approximations 
are shown in Fig. 1. The leading approximation from Eq. (31), labeled G2, deviates 
from the critical correlator strongly. Including the next-to- leading term g^ix) yields 
the approximation labeled G3, i.e. Eq. (25). A square indicates that G3 deviates 
from the critical correlator by less than 2% for t/r > 5 • 10^. If that criterion is 
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relaxed to 5%, G3 obeys it for t > IO^t. The approximation by G3 provides a first 
reasonable approximation to (t) . Including further terms of the expansion improves 
the approximation as is shown for G5 and Gy. One recognizes that proceeding from 
G5 to Gy still improves the range of applicability by one order of magnitude in time. 
We conclude that the asymptotic expansion explains quantitatively the critical decay 
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Figure 2. Time scale to in units of r for the approximation of the critical decay 
at the As-singularity of the model studied in Fig. 1 by including n orders of the 
asymptotic expansion, Eq. (31). Time scales obtained by matching G7{t) at large 
time, 35 < logj^g t < 38 are shown by crosses (x). The time resulting from 
matching the solutions at t = 10® is shown by filled circles (•). The diamonds 
(o) show the time scale resulting from matching where 0(t) =2/3. The inset 
shows to on logarithmic scale. The lines are guide to the eye. 

Matching a time scale to a-t t/r = lO^'' and using six terms of the expansion in 
Eq. (31) is not a promising perspective for fitting data. However, the expansion leads 
to a reasonable approximation also for short times. Therefore, wc may depart from the 
procedure to match to at large times and try to fit to for shorter times. Figure 2 shows 
as crosses the values obtained for to when matching the approximations at the large 
times mentioned above. We will consider two procedures for fitting. The first shall 
define a scale t'g by matching the critical correlator by the approximation at t = 10^. 
The second time scale <o is obtained from matching at 50% of the decay, i.e. for the 
time t* where 4)°{t*) = 2/3. We infer from the inset of Fig. 2 that all methods to fix 
to based on the term G2{x) alone arc off by orders of magnitude. The approximation 
G'i{x) yields the correct order of magnitude for t^ in all three approaches. Starting 
with n = 5, the scales <o and t'^ can no longer be distinguished. Therefore, matching 
the approximation at 10^ is comparable to matching a true asymptotic limit. The 
value ^0 is a better approximation for t^ than tp. 
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4. Critical correlators for one-component models at an j44-singularity 

Within tiie tiieory of tiic logaritiimic decay as presented in Rcf. [24], it is possible 
to specialize to the yl4-singularity by simply setting /is = in the final formulas. 
Different from that, the critical decay for the ^4-singularity does not follow from the 
solution for the ^a-singularity but requires a different asymptotic expansion. This 
can be inferred from the fact that all the coefficients Cmj in Eq. (26) contain fi^ in 
the denominator. However, the tricks used for finding a solution in terms of slowly 
varying functions are the same for the A4 as explained above for the A^. 



4-1. The leading contribution 

Using Eq. (106) for an A4-singularity, Eqs. (IGa) and (18) can be regrouped as 
0= 1^2(2) - /i4 5[04(t)](z) 

+ kV3(^) - M5 5[05(i)](2) 

+ n'M^) - i^&s[4>''{t)][z) 



With the Ansatz 4){t) = gm{x) = Cm/x"^ , one arrives for the terms of the first line at 
i^2{z) = C(2)(™c™/y™+i)2 + 0(l/y2™+3) .^ndS[^\t)]{z) = (c™/y™)4 + 0(l/y4™+i). 
For m = 1, the first line in Eq. (32) is of leading order 0{l/y'^) with the equation for 
the coeflacient C(2)cf = /i4cf . This results in the leading-order solution [7], 



gi{x) = ci/x, ci = VC(2)/M4 • (33) 



^.2. The leading correction 

The corrections may be rephrased in terms of a differential operator and the solution 
is straight forward as before. Since later on, only the first correction will be needed 
explicitly, it will be calculated here by the linear differential equation for the Ansatz 

m = [<t>*{t) - n /(I - n = 9i{x) + 

2y'g'{y) + Ay^g = 4 VC(2)7a^C(3)/C(2) + 3C(2)ac//.4 - M5C(2)^I ■ (34) 
This is solved in leading order by g2{x): 

g2{x) = C2 Hx)/x' , C2 = 2VC(2)/M4C(3)/C(2) + 3C(2)k/(2/.4) - M5C(2)/(2AiD .(35) 
Higher-order contributions for m ^ 3 can be written in the form 

rn—l 

= I] c„uln(a;)/x"' (36) 

1=0 

with the appropriate choice of the parameters Cm.i- Hence, the general solution for 
the critical decay at an A4-singularity in the one-component case is represented up to 
errors of order 0(ln-^"+^^ (i) as 

n 

m* =/* + (!- /*) G„(i) , G„(t) = 9ra{Ht/t^)) . (37) 

m— 1 

Because the leading order result gi{x) is of order C(l/ Int) each higher order solution 
requires the inclusion of an additional line in Eq. (32). This adds new parameters like 
/xg and k' in each step, whereas for the Aa-singularity, Eq. (26), additional parameters 
occur only in every second step of the expansion. 
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4-3. Discussion 

The results for the A4-singularity shaU be demonstrated for the kernel [6], 

m(i) = ui0(t) + i;20'W + W30^ W , (38) 

substituted mto the equatfon of motion (30) used with t = 1. The model has an 
A4-singularity at V* = (1,1,1) with /* = and coefficients fii, I ^ 4 and k being 
unity. 
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Figure 3. Critical decay <j)*{t) at the A4-singularity of the model defined 
by Eqs. (30) and (38), and the unit of time chosen such that t = 1. The 
approximations by Eq. (37) with to = 0.055 matched for G4 are labeled 
accordingly. The square and the circle mark the time where the approximation by 
G4 deviates from the solution by 5% and 10%, respectively. The triangle refers 
to a 5% deviation of G2 from the solution. The inset displays the inverse of 
[(/)*(t) — /*] and its respective approximations. 



Using up to four terms in the expansion (37), the time scale is fixed at to = 0.055. 
Successive approximations to the numerical solution are shown in Fig. 3. Again, 
the leading approximation Gi does not describe the solution. The inset shows 
[4>*{t) — /*]~^, where a decay proportional to 1/lnt would be seen as a straight line. 
Gi yields such a straight line by definition; but it has the wrong slope compared to 
the solution. The latter exhibits a straight line for t > 10^. Including the leading 
correction in G2 can account for the slope of the long-time solution. Further terms 
in the asymptotic expansion enhance the accuracy of the approximation. G4 fulfills 
the 5% criterion at t = 3 lO'^, and is in accord with the solution on the 10% level for 
t > 230. G2 intersects (j)*{t) for shorter times but deviates first from the solution by 
5% at i = 9 • 10^2. 
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5. Asymptotic expansion of the critical correlators at an As-singularity 

For the study of the general models, we go back to Eqs. (4-7). The solvability condition 
for Eq. (6) reads 

^a;j,(z) = 0, (39a) 

and the general solution can be written as 

4>git) = ag^it) + 4>q{t) . (39&) 

The splitting of (t)q{t) in two terms is unique if one imposes the convention 
'^qU* (l)q{t) = (fiit). Then, the part 0q(t) can be expressed by means of the reduced 
resolvent Rgk of aI^j^": 

S[4>q{t)]{z) = RqkJkiz) . (39c) 

The matrix Rgk can be evaluated from matrix A^]^''^ and the vectors a^., ak [35]. Let 
us emphasize that Eqs. (39a-c) together with the definitions in Eqs. (4) and (7) are 
an exact reformulation of the equation of motion (3) for states at glass-transition 
singularities. It is the aim of following calculations to express (j)q{t) recursively in 
terms of cj){t) and to show that (j){t) has the asymptotic expansion discussed in Sec. 3 
for the one-component models. The starting point is the observation that (j)q{t) is 
small of higher order than (/)(t). This is obvious, since Eqs. (7) and (39c) imply 
^q{z) = O(02) 0{$4)q) + 0{4)l). Therefore, one gets 

Jqiz) = Oi$^), (40a) 

4>q{t) ^ 0{^') . (40&) 

We assume that </> can be expanded in terms of functions gm{x) as defined in 
Eqs. (lla,fe), and show the legitimacy of this Ansatz by the success of the following 
constructions. 

5.1. Expansion up to next-to-leading order 

Substituting the splitting (396) into the inhomogeneity Jq (z) from Eq. (7) yields 

M^) = A<fl\^_ak,ak,Sm']iz) - alsm?iz) + ^(0') • (41) 

The hmction ip2{z) in Eq. (17) is of order 0{(t)^) because of Eq. (15). Therefore, 

M^) = - al)Sm-']iz) + 0{$^) . (42) 

Remembering Eq. (9) and the condition /i2 = 0, one notices that the solvability 

condition (39a) is fulfilled to order 0(0^). Hence, Eq. (39c) yields 

4>q{t) ^ Xq4;\t) + Oi$^) (43) 



'jqyuj — yvq 

with the abbreviation [24] 
A„ = R„k 



(44) 



The first step in the derivation of g-dependent corrections results in the extension of 
Eq. (396): 

0,(t) = Oqm + ^9<^'(i) + 0^(0 , (45a) 
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where 



4>'Jt)^0{$'). 



{45b) 



The next step is started by substituting the rcsuh (45a) into Eq. (7) for Jq{z). 
Terms of order 0{4P') vanish ahogether as demonstrated above and only a^V-'2(z) and 

additional terms of order 0{(jT') are left from Jq'^\z). Equation (17) is used to reduce 
products of iS-transforms to 5-transforms of products. The inhomogeneity assumes 
the form 



A 



(3)c 



Let us introduce k = 2^ and fi^ in agreement with Ref . [24] : 



(«: 



C 



[aqXq+al/2] , 



M3 = 2C - 



J2< [ 



A 



(3)c 

qki ^2^3 



qki k2 



2aqXq] 



(46) 
(47a) 
(476) 



(48) 

g3{x)-\-0{l/x^) 



Then, the solvability condition (39a) reads 

This equation was discussed in Sec. 3. The result is (f>(t) ~ g2{x) 
with the functions (72(2;) and gsix) specified in Eqs. (19) and (24), respectively. From 
Eq. (43) one infers, that 4iq{t) = 0{l/x'^). For the solution up to next-to-leading 
order, only the first term on the right-hand side of Eq. (45a) matters. However, the 
discussion of the solvability condition including the Xq(lP{t)-teTm was necessary in 
order to fix the important number ^3, which enters Eq. (48) and thereby the cited 
formulas for 52(2;) and 53 (x). 



5.2. Higher-order expansions 

After substitution of Eq. (45a) into Eq. (7) in order to extend the expansion of Jq{-. 



one can use Eq. (39c) to determine (f'qi't) up to errors of order ©(^ 
a new amplitude Yq as 



There appears 



Y„^R 



■qk 



A 



(3)c 3 



+ 2[4fe'fc2"fci^fc2 - auXk] + tizal] . (49) 



To get the last term in the curly bracket, Eq. (48) was used to express the frequency 
dependence of Jq{z) in Eq. (46) solely by S[4>{t)^]{z) . After this second reduction step, 
one gets the extension of Eq. (45 a): 



iq{t) = a,0(t) + Xq4>\t) + Yq^^'it) + (t) , 



(50 a) 



where 



<(t) = O(04). (50&) 

Here, the contribution proportional to Yq has g| as the lowest-order term, and therefore 
it is of higher order than g^. However, the calculation of the amplitude Yq is a 
prerequisite to determine the parameter /X4, which will be needed below. 

To continue, we substitute Eq. (50a) into the solvability condition (39a). The 
same tricks as before are required to yield a definition of which is consistent with the 
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equations for the one-component case. Before adding new terms from the expansion 
of Jq{z) in Eq. (7), the remaining terms of order 0{4>^) in Eq. (46) shall be collected 
from the lines with n ^ 3. A new parameter is introduced to shorten notation, 

k^2^a*aqXg, (51) 
q 

and the contribution to Jq{z) so far is k'03(z) — kS[(j)\')p2{z)- Equation (48) can be 
used to eliminate iIj2{z). With the assistance of Eq. (17), this contribution is reduced 
to Kip3{z) — ii3KS[(j)*] -\-0{(j)^). Next, the term from Eq. (7) for n = 4 is added and the 
term with k is absorbed in the definition of /i4. Then, the solvability condition reads 

= n^z) - ImS[4>^] + O(05) ^ (52) 

where the definition for the remaining parameter fi^ is 

q 

+ \X^q - A'fllXkM + 2[a,Y, - A%\\^ak,Yk,]} + . (53) 

After having defined all the necessary parameters, we see that the solution from 
Sec. 3.3 for (j){t) is consistent with the solution of the q-dependent case as formulated 
in Eq. (50a). Keeping only terms up to errors of order (1/lnt)^, one arrives at the 
asymptotic formula for the critical correlator at an As-singularity, 

</.°(t)= f° + hl{g2ix) + gsix) 

+ Mx) + K°gUx)] + [55 (x) + 2K°g2{x)g3{x)]} , (54a) 

with 

h°^{l- f°)ag , K° = X,/aq . (54&) 

The first line of Eq. (54a) expresses the factorization theorem: (/)°(i) — /° is a product 
of a first factor h°, which is independent of time, and a second factor [32(2;) + gsix)], 
which is independent of the correlator index q. Factorization is first violated in order 
1/ln^t, and only the terms with the amplitudes K° are responsible for that. The 
expansion for 0° (t) can be carried out up to order 1 / In^ t if /i4 is known. The next 
order includes ge{x) and requires knowledge of the additional parameter ^^5. 



5.3. Discussion 

As simple example for the demonstration of the preceding results, an A/ = 2 model 
shall be considered. The MCT equations for Brownian dynamics read for q = 1, 2 

Tqdt4)q{t)+4)q{t)+ j mg{t - t')dt' (pq{t' )dt' ^ , (55a) 

Jo 

mi{t) ^ vicplit) + V2(l)l{t) , (55&) 

m2{t) = V3Mt)<p2{t) . (55c) 

This is a schematic model for a symmetric molten salt [36]. The model has three 
control parameters, V = (wi,W2,W3)- The glass-transition singularities in this system 
can be evaluated analytically. There is an A4-singularity at « 24.78, and A3- 
singularities occur for > • To allow for a comparison with previous work [24] , we 
set Ti = r2 = 1 and choose the Aa-singularity for =45. 
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Figure 4. Critical decay at the As-singularity in the two-component model 
defined by Eqs. (55a-c) for = 45 and ti = T2 = 1. The rcscaled solutions 

<l)g{t) = [4>q{t) — fq]/f^qt Q = li2, arc shown for as full lines. The asymptotic 
approximation (54a) is shown dashed for g = 1 and dotted for q = 2. The 
points where the approximation deviates by 5% from the solution for q = 1,2, 
and the point where the solutions dilfer by 5% from each other are marked by 
a square, a triangle and a circle, respectively. The chain curve with label G2 
shows the leading contribution from Eq. (54a). The inset shows as full lines the 
rectification, 4>qit)~^^^ for q = 1 (lower full line) and q = 2 (upper full line). The 
q-independent part G5 of the approximation in Eq. (54a) is given by the dashed 
line. The dotted and the chain line show the leading and next-to-leading order 
approximations G2 and G3, Eq. (54a). The time scale to is 4.07 • 10"'^. 



Let US use the rescaled correlators (f>q{t) — [^^(t) — fq]/h° for the following 
considerations. The result in Eq. (54a) assumes the form <j>°{t) = G5(x) + K°Gci{x), 
with G5{x) from Eq. (31) and G^ix) = 5|(x) -f- 2g2{x)g^[x). Since G^{x) is of 
higher order than G^{x)^ Eq. (19), correlators for different q approach each other 
for sufficiently large time as is demonstrated in Fig. 4. The time t « 2 • 10^, where 
(j)2 deviates by 5% from 0°, is marked by a circle. The amplitude K° introduces the 
g-dependent corrections which are smaller for q ^ 1 than for q = 2. To evaluate 
6*5 (x) and G^lx), we determined the following parameters, I-I3 = 0.772, k ~ 0.888, 
and fi4 = 1.38. Notice, that ^3 is more than twice as for the model studied in 
Fig. 1. Since the coefficients Cm,i in Eq. (26) contain powers of yUs in the denominator, 
corrections are smaller if ^3 is larger, cf. Eqs. (27ffl)-(27c) and (28a)-(28rf). Because 
of the smaller corrections, the time scale can be matched with G5{x) between t = 10^" 
and 10^^ which is significantly earlier than for the model studied in Fig. 1. We get 
to = 4.07-10-3. 

The asymptotic approximation (54a) is shown as a dashed line for q = I in Fig. 4, 
it deviates by more than 5% from the solution if < < 10^ (□). The approximation 
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for g = 2 (dotted) deviates by more than 5% for t < 6 ■ 10^ (A). This difFercncc in 
the range of validity can be understood quahtatively by considering the g-dependcnt 
corrections of higher order in Eq. (50a), Kq[g^{x) -\-2g2{x)g4{x)] + Yqg2{x)/aq with Yq 
from Eq. (49). Both Kg and Yq/ug are smahcr for the first correlator, Yi/ai = —0.1928 
and 12/12 = 5.761, and introduce less deviations from the q-indepcndent part Ge{x) 
of the approximation in higher order. 

The (j-independent function G^^x) would lie on top of the dashed line in Fig. 4 
and is therefore shown only in the inset which also displays the critical correlators and 
the g-independent functions G2(x) and G3{x), Eq. (31). Plotting (^°(t)~^/-^ we can 
identify 1/ In^ t-behavior as straight line. The critical correlators exhibit a straight 
line starting from t w 10^. The leading approximation 6*2(0;) is a straight line as well 
but has a slope slightly larger than the solution. The first correction G3{x) resembles 
the slope of the solution but is offset from the solution by a shift of the time scale. 
This was observed before in Fig. 1. Since G'5(.t) + K°Gr,{x) was used to match the 
time scale to and as G^lx) decays faster than the g-indcpendcnt part, 65(2;) coincides 
with the solution for larger times. 

As a second example, the asymptotic laws shall be considered for the square- 
well system (SWS). This is the microscopic model for a colloid explained in Sec. 1. 
The microscopic version of MCT for colloids is used with the wave-vector moduli 
discretized to a set of M = 500 values. The structure factors that define the mode- 
coupling functional J-q in Eq. (1) are calculated in the mean-spherical approximation. 
We shall consider the same As-singularity for 6° = 0.03 as considered in previous 
studies [26,27]. The reader is referred to these papers for further details and for 
an extensive discussion of the relaxation near the specified As-singularity. For the 
evaluation of the approximation (54a), we need the correction amplitudes K° which 
arc shown in Fig. 8 of Ref. [26] and the parameters characteristic for the Aa-singularity 
under discussion, 

At3 = 0.109, K = 0.314, ^4 = 0.204. (56) 
The asymptotic approximation reads 
^°q{t) = 60A/x^ -\~ 264.71nx/x^ 

+ [3374.9 - 580.2 In a; + 870.4 W x]/x^ 

+ [-11745.7 - 27952.1 ln.T - 4452.2 In^ x + 2544.1 In^ x]/x^ 

+Kq {3643.9/a;'' + 31953.7 Inx/x^} + 0{x-^) . (57) 

The first line represents g2{x) and gsix), Eqs. (19) and (24). The second and third line 
exhibit the contributions up to 54(2;) and g^ix), Eqs. (26-28(i), which arc independent 
of the wave vector. The q-dcpendcnt correction terms appear with the prcfactor K° 
in the curly brackets; they are positive for t/to > 2.5 and monotonically decreasing 
for t/to > 5.1. 

Figure 5 shows the rescaled functions 4>°{t) for three representative wave numbers. 
At the peak of the structure factor, qd = 7, the amplitude is negative, for qd = 57.4 
the correction amplitude is close to zero, and for the wave vector qd = 172.2 the 
amplitude is positive. The functions (full lines) deviate strongly from each other 
in the window of time presented, demonstrating severe violation of the factorization 
property. If the deviations among the correlation functions for different wave vectors 
cannot be assigned to the q-dependent corrections in Eq. (57) within an accessible 
window in time, we cannot expect that Eq. (57) will be sufficient to describe the critical 
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Figure 5. Critical decay at the Aa-singularity of the square-well system (SWS) 
for the relative attraction-shell width 5° = 0.03. Pull lines show the rescaled 
correlation functions <j>g(t) = [4>q{t) — fq]/h° at V°, for the wave-vector values 
qd = 7, 57.4, and 172.2 as indicated. The unit of time is chosen so that 
1/Do = 160, with Do denoting the single particle diffusivity [26,27]. The 
dashed lines exhibit the asymptotic approximation of Eq. (57) with a time scale 
to = 4 • IQ-'^ matched in the interval t = 10*°. . . 10"^. For qd = 7.0, 57.4, and 
qd = 172.2, the correction amplitudes are K° = -1.704, -0.00224, and 4.814, 
respectively. The filled diamonds for t = 10* and t = 10^^ mark the values for 
4>q{t) for the three g-values. The inset shows <f>°{t)~^^^ for the g-values above 
from top to bottom and the g-independent approximations defined in Eq. (31) in 
the same representation, G2(t)~^/^, G3(t)^^/^ and G5{t)^^^^ , respectively. 

decay. Suppose, the critical correlators for different wave vectors are approximated 
by Eq. (57). Then, for arbitrarily chosen wave vectors qi and (72, the difference 
A[qi , (72] (t) ~ (f'qi (t) ~ 4>q2 i^) is given in leading order by the difference in the correction 
amplitudes. K°_^ "^"^^2' ^'^'-^ terms in the curly brackets in Eq. (57). From Fig. 5 we 
infer that A[qi, (72] (i) is not yet close to zero to neglect the terms in the curly brackets. 
The values of {t) for the three chosen g- values are marked by diamonds in Fig. 5 for 

t = 10^ and 10^2. We get A[7, 57.4](105) = -0.030 and A[172.2, 57.4](105) = 0.161. 
These differences are large but they correctly reflect the ordering in the values for K° 
which increase with q. From that we conclude that the treatment of the g-dependence 
in Eq. (57) is qualitatively correct. 

If the time dependence of A[qi, q2\(t) were given exclusively by the terms in 
curly brackets in Eq. (57), then the differences among the K° would explain the 

amplitudes of the decay in A[qi, q2\(t). To quantify deviations from that case we 
introduce the ratio i^[qi , q2 , q^Kt) = A[qi,q2]{t)/ A[q2,q3]{t). For t 00 this ratio 
is = — K°^)/{K°^ — K°^). Deviations from Voo indicate that higher-order 

g-dependent corrections are present in addition to the terms in Eq. (57). For the 
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g-values used in Fig. 5 we get i^oo = {Kj — K^-j a) I i^lr a ~ -^172.2) ~ 0-354. Since 
iiTgy 4 w 0, this ratio is almost equivalent to ~Kj / K^j2 2- The ratio at time t ~ 10^ 
is i^i?, 57.4, 172.2] (10^) = 0.187 and therefore deviates by 90% from v^o- Hence, we 
cannot expect Eq. (57) to describe the critical decay in Fig. 5 at that time. Att = 10^^, 
the ratio has decayed to i^[7, 57.4, 172.2] (10^^) = 0.280 which which deviates from 
Uoo by 20%. Here, the g-dependent corrections are also in reasonable quantitative 
agreement with the approximation in Eq. (57). To determine to, we use extremely 
large times. The inset of Fig. 5 displays the rescaled correlators as (^°(t)~^/^. In this 
representation, the leading term 52(2;) in Eq. (5) yields a straight line. We see that 
for large times the correlators for different q indeed come closer together and the ratio 
at t = 10"° is i/[7, 57.4, 172.2] (10^^°) = 0.341, which deviates by 4% from ly^. For the 
determination of to we use Eq. (5) for q = 7, 57 A , and 172.2 and match the asymptotic 
approximation to the numerical solutions in the interval from t = 10'"' to t = 10*^. 
This results in a value to = 4 • 10~^. For times larger than t w 10^° the numerical 
solution does no longer follow the approximation. In that region inaccuracies in the 
control-parameter values lead to deviations from the asymptotic behavior. These 
inaccuracies prevent us also from fixing more than just one digit of to- The dashed 
line in the inset labeled G5 shows the result for neglecting the last line of Eq. (57). 
This also describes the correlator for q = 57.4 where Kg is close to zero. Taking into 
account only the first line of Eq. (57) yields the dotted curve labeled G3. This curve 
is clearly inferior to G5 , but it captures the slope of the solution still better than G2 ■ 
In the large panel of Fig. 5, one can compare the critical correlators with the 
approximation by Eq. (57). For times of interest for experimental studies, the 
description is reasonable qualitatively. Especially the leading g-dependent corrections 
describe the variations seen in the correlators down to relatively short times. The 
accuracy of the approximation that was demonstrated for the schematic models in 
Figs. 4 and 1 is far better than seen in Fig. 5 for the SWS. This difference is mainly due 
to different values of the parameter fi^ that characterizes the various Aa-singularities. 
For the two-component model we had fi^ = 0.77 and for the one-component model 
there was = 1/3. The small value /is = 0.109 for the SWS implies slow convergence 
of the asymptotic expansion. Therefore, a quantitative description by Eq. (57) is 
possible only for times exceeding considerably the ones shown in Fig. 5. 

6. Asymptotic expansion of the critical correlators at an A4-singularity 

6.1. Expansion up to next-to-leading order 

The calculation of the critical correlator at the A4-singularity is so involved, that we 
restrict ourselves to the leading and next-to-leading order result. The Eqs. (40a-53) 
remain valid, and Eqs. (49) and (53) simplify because /is = 0. The difficulty comes 
about because /is, which enters Eq. (35), has to be determined. This requires the 
extension of Eq. (50a), and thereby there appears a further amplitude. The additional 
amplitude Zq is obtained by also including terms with n ^ A from Eq. (7). Applying 
the same manipulations as above, one arrives at <t>'q{t) ~ Zqcf)^ + 0{4/') with the 
amplitude 

Zq = Rqk{[Af^'^k^k^j.^ak^ak.,ak,ak^ - al] + 3[A{,^fc^\^j.^afciafe2 - a^Afe] 

+ [^ifcl fc. A,., - A^] + 2 au, Yu. - afcFfe] + y^^al } . (58) 
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Introducing the third g-dependent correction, the solution assumes the form 

0,(0 = ag${t) + Xq^'^it) + Yq$^{t) + Zq^*{t) + 0{^^) . (59) 

Collecting all terms of order 0(0^) after including also the line n = 5 from Eq. (7), 
one gets from the solvability condition (39 a): 

q 

+ 4[afe^fe - ^tk^k2k3kiO-ki"'k2ak3Xki] 

+ S[akXl + alYk - A'j^^^j^^^^iak.Xk^Xk, + ak.ak^Yk:,)] 

+ 2[XkYk + akZk - A'^;^l';^^{Xk,Yk2 + ak.Zk,)]} + . (60) 

Summarizing, the asymptotic solution for the critical decay at an ^4-singularity 
in next-to-leading order reads 

<P*,{t) = /; + K {51 (^) + [92{x) + Klgl{x)]] . (61) 

Here, in analogy to Eq. (546), the critical amplitude is h* = {1 — f*)aq and the 
correction amplitude is given by K* = Xq/uq. The factorization theorem is obeyed 
by the leading-order term only. Contrary to what was found in Eq. (54a) for the 
behavior at the ^s-singularity, already the leading correction term g2 is modified by 
the g-dependent term K*gf(x) of the same order. The higher-order contributions 
enter the curly brackets in Eq. (61) as 53(2;) -f 2gi{x)g2{x)Xq/aq + gi{x)Yq/aq and 
34(2;) -f gl{x)Xq/aq -I- 2gi[x)g'i{x)Xq/aq -I- 'igl{x)g2{x)Yq/ aq + g\{x)Zq/aq. However, 
33 (x) requires the evaluation of the parameters /xg and k', .94 (.t) needs and k!' . 

6.2. Discussion 

Figure 6 shows the critical decay at the A4-singularity of the two-component model 
defined in Eqs. (55a-c). The parameters for the evaluation of gi{x) and 32(2;) are 
/Z4 = 1.53, ^5 = 0.962, and n = 0.386. We use again the rescaled correlator 
4'*q{t) = Wq{t)^ fq]lh*q and chcck first the validity of the factorization in Eq. (61) in the 
form (j)*q{t) = G2{x) + KqG2{x) where 6*2(2;) = gi{x) -I- 32(2:) and G2{x) = 3i (2:). The 
time, where the solutions for g = 1, 2, differ by 5% is only reached aX t k, 10^^ . The 
circle marks the point where the deviation is still 10% at t = 10^^. We can then use 
the approximation (61) to fix the time scale to to = 2.0 which then yields the dashed 
and dotted curves for g — 1, 2, accordingly. For q ~ 1 this approximation deviates by 
5% from the solution at t w 8.2 ■ lO"* (□). For q = 2 we find t w 1.8 ■ 10*^ (A). This 
is plausible when appealing to the g-dependent higher-order correction in Eq. (59), 
which incorporates in addition to drastically different values for Kq also the values 
Fi/oi = —0.579 and 12/02 = 3.76. A rectified representation of the critical decay 
and the approximation in the inset shows again the leading-order Gi {x) (dotted) as a 
straight line of different slope than the solution (full lines) and the second correction 
G2{x) (dashed). In this plot, the critical correlators for different q are still significantly 
different in the entire window. But Eq. (61) can account for that difference as is shown 
by the good agreement of the curve labeled G2+G2K2. The latter describes the second 
correlator where the deviations due to the correction amplitudes are largest. 

We now turn to the yl4-singularity of the SWS. For the application of Eq. (61) 
we need the parameters characterizing the j44-singularity, 

/X4 = 0.131, K = 0.243, /Lis = 1.21. (62) 




Figure 6. Rescaled critical decay <j)* = [<j>g{t) — fq]/h* at the A4-singularity in 
the two-component model defined in Eqs. (55a-c) (full lines). The asymptotic 
approximations, Eq. (61), for g = 1,2, are represented by the dashed and dotted 
curve, respectively. For q = 1 (□) and q = 2 (A), the points are marked where the 
solution and the approximation deviate by 5%. An additional point is indicated 
where the solution for q = 2 differs from the one for g = 1 by 10% (O)- The 
inset displays the rectified representation of the solutions for q = 1 (lower full 
line) and q = 2 (upper full line) together with the g-independent parts of the 
approximations, Gi and G2, cf. Eq. (37), and G2 + K2G2 (see text). The time 
scale to = 2.0 was matched for t = 10^° . . . 10^^. 



The rather small value of ^l^ generates particularly large coefficients in the expansion 
of the critical decay in Eq. (36) where appears in the denominators. This feature 
is quite the same as mentioned above for the Aa-singularities. The asymptotic 
approximation in Eq. (61) yields for the critical decay of the rescaled correlators: 

0;(i) = 2,M/x-bQ.7\nx/x^ + l2.bKl/x^ + 0{x-^) . (63) 
We chose again values for q where K* is negative, almost zero and positive. Figure 7 

demonstrates that the factorization is strongly violated. Comparing the solutions (/)*(t) 
for t = 10^ we find a ratio defined as in the previous section of z^[7, 32.2, 39.8] (10^) = 
1.439 which is more than 30% off the ratio for the correction amplitudes 1^00 = 2.185. 
At t = 10^2 zy[7, 32.2, 39.8] (10^2) = 1.723 which achieves 20% accuracy. So 

the critical decay at the A4-singularity shown in Fig. 7 is in qualitative accord with 
Eq. (63) with respect to the variation in q. However, due to the small value of /X4, 
the differences among the correlators for different q do not decay fast enough to allow 
for a consistent determination of to for the maximum value in time that could be 
reached. Numerically we find i/[7, 32.2, 39.8] (lO^^S) ^ 2.076 which is stiU 5% off from 
z/qo, and 4)*q{t) itself deviates from zero by 5%. This illustrates drastically the enormous 
stretching at the A4-singularity. 

The inset of Fig. 7 demonstrates that the critical decay (/)J(t) is qualitatively 




Figure 7. Critical decay at the j44-singularity of the SWS for qd = 7.0 (dashed), 
32.2 (full line), 39.8 (dotted). The correction amplitudes are Kg = -1.81, -0.04, 
and 0.77, respectively. The filled diamonds mark the values at i = 10^ and 
t = 10^^ where the ratios u{t) are 1.44 and 1.72, respectively (see text). The inset 
replots the curves from the full panel in the same linestyle and shows the first 
term of Eq. (63) labeled Gi and the law ln(t/r)~^/'^ labeled A^, both with an 
arbitrary time scale. 



different from the leading order 1/ Int-law for t ^ 10^". For the As-singularity in Fig. 5 
it was still possible to argue that curve G2 is in accord with the decay qualitatively 
at least for large times and to attribute deviations for shorter times to the proximity 
of the A4-singularity. Figure 7 does not allow for such an interpretation. The curves 
l/(f>*{t) have a slope smaller than 1/Gi over the complete window in time and imply a 
slower decay than given by the leading order in Eq. (63). If fi^ was zero, the singularity 
would be of type A5. The leading order critical decay at such butterfly singularity is 
ln(i)~^/^. This law is added in the inset as chain line labeled A^. Indeed, it explains 
the data qualitatively. Hence, the shortcomings of the asymptotic expansion at the 
A4-singularity in the SWS result from the small value of fi4. 

To check if the value for ^4 is exceptionally small for the SWS, the calculation 
was repeated for the hard-core Yukawa system as introduced in Ref. [16]. We find the 
even smaller value 114 = 0.080. Therefore, the small value of /i4 seems to be typical 
for systems with short-ranged attraction. 

7. Summary 

The asymptotic expansion for large times of the critical decay of correlation functions 
at higher-order glass-transition singularities has been elaborated. These decays 
can be considered as the analogue of the t~''-law expansion for the correlators at 
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the liquid-glass transition. The latter as well as the higher-order singularities are 
obtained as bifurcations of type Ai, I ^ 2. The A; -singularity and especially the 
critical decay law at the singularity is characterized by a number /x;. For the A2- 
singularity of the liquid-glass transition, this characteristic number determines the 
so-called exponent parameter A = 1 — /i2, which specifies the critical exponent a via 
A = r(l — a)^/r(l — 2a). For /i2 = or A = 1, one gets a = and the asymptotic 
expansion in terms of powers becomes invalid. A higher-order singularity An is 
encountered, defined by > while /x; = for Z < n. 

For an Aa-singularity, the critical decay is given by an expansion in inverse powers 
of the logarithm of the time, starting with 1/ In^ t. The convergence of the asymptotic 
expansion is the better the larger is fi^. The result for the general models in Eqs. (54a) 
and (26) adds probing-variable dependent correction terms to the one-component 
result. These can be expressed by terms from the one-component solution and 
correction amplitudes. The leading correction amplitude Kq is the same function of 
the MCT-coupling constants as found earlier for the logarithmic decay-law expansions 
[24]. Since the vertex is a smooth function of the control parameters, these correction 
amplitudes are smooth functions as well. Therefore, also for the general case, the 
range of validity for the asymptotic expansion is determined by the characteristic 
parameter fi^. If fi^ is small, the quality of the fit by the asymptotic expansion is less 
satisfactory than for larger ^3. Generically, larger ^3 can be obtained by extending the 
corresponding glass-glass-transition line deeper into the glassy region and hence having 
the As-singularity further separated from the liquid regime. Thus, the dynamics 
influenced by an As-singularity seen in the liquid regime is either connected to a rather 
small fX3, or it is strongly influenced by a crossing of different liquid-glass-transition 
Hues [27]. 

For /i3 =0, an A4-singularity is found; the expansion for one-component models 
in Eqs. (26-28d,31) becomes invalid and has to be replaced by Eqs. (36) and (37). 
The general solution in Eq, (61) has similar properties as mentioned above for the A3- 
singularity. Now it is the characteristic parameter /i4 that determines how satisfactory 
the approximation can be. While ^4 = 1 in Fig. 3 and ^4 = 1.53 in Fig. 6 allows for a 
description in the schematic models considered, the small parameter fj.4 « 0.1 in the 
microscopic models for systems with short-ranged attraction prevents the application 
of the asymptotic formula. 

An understanding of the critical decay law is a prerequisite for estimating the 
range of validity of the Vogel-Fulcher-type laws which describe the asymptotic limit of 
the time scale of the logarithmic decay laws near the higher-order singularities [7] . For 
the two-component model analyzed above, the asymptotic limits were demonstrated 
for reasonable windows in time [25]. For the mentioned colloid models, the small 
values of the characteristic parameters and fi^ together with the manifest violation 
of the factorization property restrict such laws to unreasonably long times. 
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